Decoherence as a relative phenomenon: a generalization of the spin-bath model 
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In this paper we propose a closed-system perspective to study decoherence. ^From this perspective 
we analyze the spin-bath model as presented in the literature, and a natural generalization of that 
model. On the basis of the results obtained from that analysis, we argue that decoherence may be 
understood as a phenomenon relative to the partition of a closed system, selected in each particular 
case. This viewpoint frees the decoherence program from certain conceptual difficulties derived from 
its open-system perspective. We also argue that the usual picture of decoherence in terms of energy 
■ dissipation is misguided. 
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' In the literature on quantum mechanics, 'decoherence' refers to the quantum process that turns a coherent pure 
state into a decohered mixed state, which is diagonal in a well defined basis. The phenomenon of decoherence is 
essential in the account of the emergence of classicality from quantum behavior, since it explains how interference 
vanishes in an extremely short decoherence time. 

The orthodox explanation of the phenomenon is given by the so-called environment-induced decoherence (EID) 
approach ([![, 0], 0]), according to which decoherence is a process resulting from the interaction of an open 
quantum system and its environment. In particular, by studying different physical models, it is proved that the 
reduced state pg(t) = TtePse^) °f the open system rapidly diagonalizes in a well defined pointer basis, defined 
case by case but not in general, which identifies the candidates for classical states. The EID approach has been 
extensively applied to many areas of physics, such as atomic physics, quantum optics and condensed matter. In 
particular, the study of decoherence has acquired a great importance in quantum computation, where the phenomenon 
l— — '■ of decoherence represents a major obstacle to the implementation of information processing hardware that takes 
advantage of superpositions. In spite of its impressive practical success, from a conceptual viewpoint the EID approach 
still faces a difficulty derived from its open-system perspective: the problem of defining the system that decoheres. 
£T) ' ^From the einselection view, the split of the Universe into the degrees of freedom that are of direct interest to 
the observer —the system— and the remaining degrees of freedom —the environment— is absolutely essential for 
decoherence. However, since the environment may be external (a "bath" of particles interacting with the system of 
interest), internal (such as collections of phonons or other internal excitations) or a combination of both cases, the 
t — ■ EID approach offers no general criterion for deciding where to place the "cut" between system and environment. In 
many cases, the lack of such a general criterion leads to the need of assuming in advance the observables that will 
behave classically. For instance, in cosmology the usual strategy consists in splitting the Universe into some degrees 
of freedom representing the "system" of interest, and the remaining degrees of freedom that are supposed to be non 
accessible and, therefore, play the role of an internal environment. In quantum field theory, when it is known that the 
background field follows a simple classical behavior, the scalar field is decomposed according to 4> = <fi c + 4> q , where 
the background field <fi c plays the role of the system and the fluctuation field (j> q plays the role of the environment 
(see Q). Zurek concedes that this absence of a general criterion to discriminate between system and environment is 
a serious difficulty of his proposal: "In particular, one issue which has been often taken for granted is looming big, as 
a foundation of the whole decoherence program. It is the question of what are the "systems" which play such a crucial 
role in all the discussions of the emergent classicality. This issue was raised earlier, but the progress to date has been 
slow at best" (see @, p. 122; for a discussion of this point, see Q). 

The main purpose of this paper is to argue that this "looming big" problem is actually a pseudo-problem, which is 
simply dissolved by the fact that the split of a closed quantum system into an open subsystem and its environment 
is just a way of selecting a particular space of relevant observables of the whole closed system. But there are many 
different spaces of relevant observables depending on the observational viewpoint adopted. Therefore, the same 
closed system can be decomposed in many different ways: each decomposition represents a decision about which 
degrees of freedom are relevant and which can be disregarded in each case. Since there is no privileged or "essential" 
decomposition, there is no need of an unequivocal criterion for deciding where to place the cut between "the" open 
system and "the" environment. On this basis, we will show that the usual picture of decoherence in terms of energy 
dissipation from the open system to the environment can no longer be sustained. Summing up, decoherence is a 
phenomenon relative to the relevant observables selected in each particular case. The only essential physical fact 
is that, among all the observational viewpoints that may be adopted to study a quantum system, some of them 
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determine subspaces of relevant observables for which the system decoheres. 

For the purpose of this argumentation, the paper is organized as follows. In Section II, by means of the concept of 
tensor product space, it will be shown how the split of a whole closed system into an open system and its environment 
can be understood as the selection of a space of relevant observables. In Sections III, IV and V, the well-known 
spin-bath model studied by the EID approach is presented from the perspective of the previous subsection, and 
physically relevant numerical simulations are obtained. In Sections VI, VII and VIII a generalization of the spin-bath 
model is presented and solved by computer simulations; this task will allow us to compare the results obtained for 
different ways of splitting the entire closed system into an open system and its environment. In Section IX, the results 
obtained in the previous sections are discussed from a conceptual viewpoint in order to argue for the relative nature 
of decoherence —which dissolves Zurek's "looming big" problem— and for the rejection of the usual description of 
decoherence in terms of energy dissipation. Finally, in Section X we draw our conclusions. 



II. SELECTING THE RELEVANT OBSERVABLES 



As it is well-known in the discussions about irreversibility, when a —classical or quantum— state evolves unitarily, 
it cannot follow an irreversible evolution. Therefore, if a non-unitary evolution is to be accounted for, a further 
element has to be added, precisely, the split of the maximal information about the system into a relevant part and an 
irrelevant part: whereas the irrelevant part is discarded, the relevant part may evolve non- unitarily. This idea can be 
rephrased in operators language. Since the maximal information about the system is given by the space O of all its 
possible observables, then we restrict that maximal information to a relevant part by selecting a subspace Or C O 
of relevant observables. The irreversible evolution is the non-unitary evolution viewed from the perspective of those 
relevant observables. 

As emphasized by Omnes (Q, @), decoherence is a particular case of irreversible process. Then, the selection 
of the subspace Or C O is always required in decoherence. In the case of the EID approach, the selection of Or 
amounts to the partition of the whole closed system U into the open system S and its environment E (see [Io|). In 
fact, let us consider the Hilbert space H of the closed system U, H = Hsd> He, where Hs is the Hilbert space of S 
and He the Hilbert space of E. The corresponding von Neumann-Liouville space off/ is C — H ®H = £s® Ce = O, 
where Cs = Hs® Hs and Ce = He® He- In the EID approach, the relevant observables are those corresponding to 
the open system S: 

Or = O s ®I e eO R cO (1) 

where Os G Cs and Ie is the identity operator in Ce- The reduced density operator p s (t) of S is defined by tracing 
over the environmental degrees of freedom, 

Ps (t) =Tr EP (t) (2) 

The EID approach studies the time-evolution of Ps(t) governed by an effective master equation; it proves that, under 
certain definite conditions, pg(t) converges to a stable state p Sae : 

Ps(t) — ► Ps* (3) 
But we also know that the expectation value of any Or G Or in the state pit) of U can be computed as 

(Or) pW = Tr (p(t)(O s ® I E )) - Tr (p s (t) O s ) = (O s ) Ps{t) (4) 

Therefore, the convergence of Pg{t) to p Slf implies the convergence of the expectation values: 

(Or) pW = (Os) Ps{t) — > (O s ) Pst = {0 R ) Pt (5) 

where p % is a final diagonal state of the closed system U, such that p St = Tte p* (see eq. ©; for details, see [101). 
This means that, although the off-diagonal terms of p(t) never vanish through the unitary evolution, decoherence 
obtains because it is a coarse-grained process: the system decoheres from the observational point of view given by any 
observable belonging to the space Or. 

When viewed from this closed-system perspective, the discrimination between system and environment turns out 
to be the selection of the relevant observables. By following papers [ll[ and [l^], we will use the expression ' 'tensor 
product structure' (TPS) to call any factorization H = Ha^Hb of a Hilbert space H, defined by the set of observables 
{Oai <8> Is, Ia ® Osi}, such that the eigenbases of the sets {O^i} and {Osi} are bases of Ha and Hb respectively. If 
H corresponds to a closed system U, the TPS H — Ha®Hb can be viewed as representing the decomposition of U into 
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two open systems Sa and Sb, corresponding to the Hilbert spaces Ha and Hb respectively. In turn, given the space 
O = H <8> H of the observables of U, such a decomposition identifies the spaces O a = Ha <8> Ha and = Hb ® 
of the observables of the open systems SU and Sb, such that Oa 8IbCO and 1,4 eg) Ob C C Once these concepts 
are considered, the selection of the space Ob. of relevant observables in the EID approach turns out to amount to the 
selection of a particular TPS, H = H S ® He, such that O r = O s ® l E C O = H ® H. 

In this paper we will consider the particular case where the closed system U is composed of n spin-1/2 particles Pj, 
each represented in its Hilbert space He 

n 

H = H!®H 2 ® ...®W„ = (g)Wi (6) 

It is quite clear that the system U can be decomposed into two subsystems S and E in different ways, depending 
on which particles are considered as the open system S. For instance, if the particle P\ is the open system S, the 
corresponding TPS reads 

H = Hs®H E = {H 1 )® f& w ^ ( 7 ) 
In turn, if the particle Pk, with 1 < k < n, is viewed as the system S, the corresponding TPS is 

H = H s ®H E = {H k )® (8) 

\i#fc / 

But we can also define the system S as composed of more than a single particle; for instance, if the particles Pj, with 
j = 1 to m < n, are the system S, the TPS in this case reads 



H = Hs®H E = \<S> n 3 I® «i (9) 

\j=l / \i=m+l / 

In the following sections we will study the phenomenon of decoherence for different partitions of the whole closed 
system U. 



III. THE SPIN-BATH MODEL 
A. Presentation of the model 

The spin-bath model is a very simple model that has been exactly solved in previous papers (see [l[). Here we will 
study it from the closed-system perspective presented in the previous section. 

Let us consider a closed system U = PUPi U. . . UP/v = PU (U-^P,), where (i) P is a spin-1/2 particle represented 
in the Hilbert space Hp, and (ii) each Pj is a spin-1/2 particle represented in its Hilbert space Hi- The Hilbert space 
of the composite system U is, then, 

h = h p ®(^$h}j (10) 

In the particle P, the two eigenstates of the spin operator S P1 ? in direction 1? are |f|*) , | -(]•): 

£p,v \t) = \ it) s Ptlt m = -\ \\) (ii) 

In each particle Pj, the two eigenstates of the corresponding spin operator S i in direction ~~v are ||j) , 

IT*) - \ %) \u) = -1 \u) (12) 
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Therefore, a pure initial state of U reads 

|Vo) = (« + b W)) ® (fg)(a 2 | T<> + Pi\ U))j (13) 

where |a| 2 + \b\ 2 = 1 and |a;| 2 + |/3J 2 = 1. If the self-Hamiltonians Hp of P and ffj of Pi are taken to be zero, 
and there is no interaction among the Pi, then the total Hamiltonian H of the composite system U is given by the 
interaction between the particle P and each particle Pi (see [l3l|): 

1 N 

i=l 

where Ij = | ) (Tjl + l-Lj) (Ijl is the identity operator on the subspace Ttj. Under the action of H, the state | V^o) 
evolves into 

m))=a\t)\£ 1[ (t)}+b^)\S^t)) (15) 

where 

N 

\£f(t)) = M-t)) = (g) (a, e~^ 2 IT,) + P % e^' 2 \U)) (16) 
t=i 



^ (It,) (Til-UO U<D® 




(14) 



B. Computing the expectation values 

The space O of the observables of the composite system U can be obtained as O — Op ® where Op is 

the space of the observables of the particle P and Oi is the space of the observables of the particle P$. Then, an 
observable O ^ O = TL ®TL can be expressed as 

N 

= P ®{(g)O l ) (17) 

i=l 

where 

Op = s n |fr) (fr| + « M |fr) m + s u |4) (1H + s M |4) (4| e Op (is) 
o, = efj | TO (Tt I + 4? I i<) & I + 4t I W (t< I + 4; I TO U< I e o< (19) 

Since the operators Op and Oi are Hermitian, the diagonal components s^, S|i4, eW, eW are real numbers, and the 

off-diagonal components are complex numbers satisfying s^j. — s^, e^? = 4*** . Then, the expectation value of the 
observable O in the state \i/>(t)) of eq. (JTSJ) can be computed as 

(0) m = (\a\ 2 s n + \b\ 2 s u )T (t) + 2Re[ab*s u r 1 (t)] (20) 

where (see [l3j ) 

N 

r (*) = II [NMt + IAI 2 4l + 2Re(a i #ege i * t )] (21) 

i=l 
iV 

rx(t) = II [klM?^ + Iftl^U 6 "^' + 2 Re(a. ( /?*£«)] (22) 

i=l 

By contrast to the usual presentations, we will study two different decompositions of the whole closed system U into 
a relevant part and its environment. 
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IV. THE SPIN-BATH MODEL: DECOMPOSITION 1 
A. Selecting the relevant observables 

In the typical situation studied by the EID approach, the open system S is the particle P, and the remaining 
particles Pi play the role of the environment E: S = P and E = u£L]Pj. Then, the TPS for this case is 



H = Hs®H E = {Up)® (® n ^J 



(23) 



Therefore, the relevant observables Or of the closed system U are those corresponding to the particle P, and they 
are obtained from eqs. (flTf . ([18]) and (fT9|) . by making = ej^ = 1 and = 0: 

On = O s ®I E = ( s *s'\s){A ) ® f® 1 *) ( 24 ) 

\s,s'=1>,^ / \*=1 / 

The expectation value of these observables in the state °f e 1- $M ^ s given by 



(0*)*(t) = l fl | 2 %tf + l^l 2 «w + 2 Re [ afe * S M K*)] ( 25 ) 



where 

N 

r{t) = (^(t)>|£fr(t)> = I] (N* ^ iffi * + Iftl 2 e^) (26) 

i=l 

and, then, 

N 

\r(t)\ 2 = Ild^l 4 + lAI 4 + Sla.H^I 2 cos2 5i i) (27) 
This means that, in eq. ([20]), r (t) = 1 and Ti(t) = r(t). 

B. Computing the behavior of the relevant expectation values 

In order to know the time-behavior of the expectation value of eq. (|25[) . we have to compute the time-behavior of 
r(t). If we take \aii\ and as random numbers in the closed interval [0, 1], such that \cti\ 2 + |/3J 2 = 1, then 

max(|a 4 | 4 + |/3 4 | 4 + 2|a l | 2 |/3 l | 2 cos2. 94 i) = ((|a 4 | 2 + I/3J 2 ) 2 ) = 1 

min(|a l | 4 + | / 3 l | 4 + 2|a l | 2 |/3 4 | 2 cos(2^)) = ((M 2 - l&| 2 ) 2 ) = (2 - l)' (28) 

Therefore, (|ai| 4 + I/3J 4 + 2\ai\ 2 \f3 i \ 2 cos2git) is a random number which, if t 7^ 0, fluctuates between 1 and 
^2 \ai\ 2 — 1^ . Let us notice that, when the environment has many particles (that is, when N — > 00), the statis- 
tical value of the cases |aj| 2 = 1, |/3J 2 = 1, | 2 = and \[3 i \ l — tends to zero. In this situation, eq. (f2"T|) for |r-(t)| 2 
is an infinite product of numbers belonging to the open interval (0, 1); as a consequence, 

lim r(t) = (29) 

N—^co 

If we know that, for N — > 00, r(t) = for any t ^ and r(0) = 1, it can be expected that, for N finite, r(t) will 
evolve in time from r(0) = 1 to a very small value. In order to obtain the time-behavior of r(t), different numerical 
simulations have been performed, where the random \oti\ were obtained from a random-number generator, and the 
were computed as = 1 — \ai\ . The value of the gt and the time-interval [0, to] for the computations were 
stipulated. The time-interval [0, to] was partitioned into intervals At — to/200, and the function |r(t)| 2 was computed 
at times tk = kAt, with k = 0, 1, 200, according to eq. (|2"T|) . 




FIG. 1: Figure 1: Evolution of \r(t)\ 2 for g t = 200Hz, and TV = 10 7 (dot line), TV = 10 8 (dash line) and TV = 10 9 (solid line), 
with t = 6.10" 6 s. 

Simulation (a): The computations were performed with TV = 10 7 , TV = 10 8 and TV = 10 9 . All the gi were taken 
to have the same value. The value gi = 400Hz was selected as a reference value on the basis of the measurement of 
the coupling constant in typical models of spin interaction ([14]). Figures 1, 2 and 3 show the time-evolution of |?"(t)| 2 
with gi — 200Hz, gi — 400Hz and g. L — 800Hz, respectively. 

In these figures we can see that, as expected, (i) for any given value of gi, decoherence is faster as TV increases, 
and (ii) decoherence is faster as gi increases, that is, as the interaction between the particle P and each particle P, is 
stronger. 

Simulation (b): Again, the computations of |?"(t)| 2 were performed with TV = 10 7 , TV = 10 8 and TV = 10 9 , but 
now the values of the gi were obtained from a random-number generator in the interval [0,800-ffz]. The results are 
shown in Figure 4 for a given distribution of the values of the gi (since the plots obtained for different distributions 
were indistinguishable) . 

Again, the figure shows that decoherence is faster as TV increases. Moreover, if we compare Figure 4 (<?,; G [0, 800Hz]) 
with Figure 2 (gi = 400Hz), we can see that the random character of the gi improves the "efficiency" of decoherence: 
the decoherence time in the case of random is shorter than in the case of constant gi. 

Simulation (c): In order to obtain a physically meaningful value of the decoherence time, we performed the 
computation for TV = 10 20 (closer to the Avogadro number) with the following strategy: 

- All the gi were taken to have the same value, under the reasonable assumption that the particles of the environment 
are all of the same nature and, therefore, all of them interact with the particle P in the same way. As explained 
above, the value of <?.; — 400Hz was selected on the basis of the typical interaction between spins. 



- The upper limit to of the time interval [0, to] was taken as 2 10 8 s, in order to show the physical dynamics of the 
phenomenon. 
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FIG. 2: Figure 2: Evolution of \r(t)\ 2 for g t = 400#z, and TV = 10 7 (dot line), TV = 10 8 (dash line) and TV = 10 9 (solid line), 
with t = 3.10" 6 s. 

- Of course, the running time of the computing process for TV = 10 20 is unattainable. Nevertheless, different runs 
of the computing process (with different random values of the \cti\ 2 and the |/3J 2 ) show that, for TV = 10 10 , 
the resulting plots are completely indistinguishable. Therefore, the result for TV = 10 20 can legitimately be 
computed by multiplying the result obtained for TV = 10 10 by itself 10 10 times. 

It is interesting to see how the decoherence time decreases as TV increases. For this purpose, we have used the 

strategy of multiplying |r(i)| 2 corresponding to TV = 10 10 by itself 10° times: (|Ki)| 2 ) 10 corresponds to the number 

of particles TV = 10 10+a . Figure 5 shows the time-evolution of (^(i)! 2 ) 10 , corresponding to TV = 10 10+a , for a = 0, 
1,2 and 3. 

The decoherence time was obtained by fitting the curve for TV = 10 20 to an exponential, and by computing the 
characteristic time of the exponential. The decoherence time so obtained was to = 10 _13 s, as empirically measured 
(see [Hj])- Let us notice in eq. (|27p that |r(t)| 2 comes back to its initial value when 2git = 27r; then, the —recurrence— 
Poincare time is tp = it/ gi. Although there is no strict final relaxation due to the discrete nature of the model, the 
relaxation time can legitimately be taken as t# < tp/2 ~ 10 -3 s. This means that, as expected, the decoherence time 
is many orders of magnitude shorter than the relaxation time. 



V. THE SPIN-BATH MODEL: DECOMPOSITION 2 
A. Selecting the relevant observables 



Although in the usual presentations of the model the open system of interest is P, we can conceive different ways 
of splitting the whole closed system U into an open system S and its environment E. For instance, we can decide to 
observe a particular particle Pj of what was previously considered the environment, and to consider the remaining 



particles as the new environment, in such a way that S 



P 3 and E 



PU(U 
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Pi). The total Hilbert space of 
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FIG. 3: Figure 3: Evolution of \r(t)\ 2 for g t = 800Hz, and TV = 10 7 (dot line), TV = 10 8 (dash line) and TV = 10 9 (solid line), 
with t = 2.10" 6 s. 



the closed composite system U is still given by eq. (fTU|) , but in this case the corresponding TPS is 





( 




H = H S ®H E = {Hi)® 


Hp® 








Km J J 



and the relevant observables Or of the closed system U are those corresponding to the particle Pj : 





(s \ 








) 




i=l J 









where (see eq. ([15])) 



o Pj = 4? I TiXTj I + 4? I ii>Ui I + 4? I iiXtj I + 4? I TjXXj I 



(30) 



(31) 



(32) 



(?) (?) (7) 

Ip is the identity operator on the subspace Hp, and the coefficients , , are now generic. The expectation 



value of the observables Or in the state \ip(t)) of eq. (fTo*)) is given by 

(o R ) m = (m\o Rj \m) = i«/4? + l^l 2 4? + 2Re ( a ^4T e ^*) 



(33) 
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FIG. 4: Figure 4: Evolution of \r{t)\ 2 for g t G [0,800Hz], and TV = 10 7 (dot line), N = 10 8 (dash line) and N = 10 9 (solid 
line), with to = 3.10" 6 s. 



B. Computing the behavior of the relevant expectation values 



In order to know the time-behavior of the expectation value of the On, we have to compute the time-behavior of 
the third term of eq. (|33p , which can be rewritten as 



2 Re 



(a j (3 



cos (g.jt + arg (a^e^JJ 



(34) 



In this case, a numerical simulation is not necessary to see that eq. (|34|) is an oscillating function which, as a 
consequence, has no limit for t — > oo. This result is not surprising, but completely reasonable from a physical point 
of view. In fact, with the exception of the particle P, the remaining particles of the environment E are uncoupled to 
each other: each Pi evolves as a free system and, for this reason, E is unable to reach a final stable state. 



VI. A GENERALIZED SPIN-BATH MODEL 
A. Presentation of the model 



Let us consider a closed system U = A U B where: 

(i) The subsystem A is composed of M spin-1/2 particles A,-, with i = 1, 2, M, each one of them represented in 
its Hilbert space H.a ■ In each Aj, the two eigenstates of the spin operator S^. in direction v are and 



S Ai ,Tt \h) = 2 s A t ,ir \h) = - j \h) 



(35) 
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FIG. 5: Figure 5: Evolution of (\r(t)\ 2 ) 10 for g t = AQQHz and a = (solid line), a = 1 (dot line), a = 2 (dash line), a = 3 
(dot-dash line), with to = 2.1CT 8 s. 

M 

The Hilbert space of A is Ha = H-Ai ■ Then, a pure initial state of A reads 

i=l 

M 

|V>a> - (g) K + &i l^» , with |a 4 | 2 + |fo 4 | 2 = 1 (36) 

i=l 

(ii) The subsystem B is composed of AT spin- 1/2 particles Bfc, with k = 1, 2, JV, each one of them represented in 
its Hilbert space Ti.B k - In each Bk, the two eigenstates of the spin operator S Bk in direction ~~v are and 
life): 



S Bfe ,vlT fe ) = ^|T fe ) S BfciV |U) = ~m) (87 



JV 

The Hilbert space of B is Hb = ® Then, a pure initial state of B reads 

fc=i 

AT 



\4>b) =0KITfc)+/3fe|lfe)), with K| 2 + |/? fe | 2 = l (38) 
fe=i 



The Hilbert space of the composite system {/ = A U B is, then, 



h = Ha®h b = (6§n A )\ ® fe)« Bfe l (39) 
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(a) 



(b) 



FIG. 6: Figure 6: 
B (white circles): 



Schema of the interactions among the particles of the open system A (grey circles) and of the open system 
(a) original spin-bath model (M — 1), and (b) generalized spin-bath model (M 7^ 1). 



Therefore, from eqs. (|36|) and (|38|) . a pure initial state of U reads 

M \ 



IV>o) = IV>a) ® \$b) 



' {<H %) + h |^» ® (g) (a* I TO + /3 fc |U» 
j / \> 1 / 



(40) 



As in the original spin-bath model, the self-Hamiltonians Ha^ and Hs k are taken to be zero. In turn, there is no 
interaction among the particles Ai nor among the particles B^. As a consequence, the total Hamiltonian H of the 
composite system U is given by (see eq. {HI) 



H = H A <E)H B 



1 









Ms 







gk (\U) (Tfcl-Uk) U*D® 




(41) 

where 1^. = (t|-j| + j-IJ-j) (JJ-j| is the identity on the subspace Ha, and I_b, = (|;| + ||;) (|/| is the identity on 
the subspace Hb, . Let us notice that the eq. (fT3)) of the original model is the particular case of eq. (|4*Tj) for M = 1. 
This Hamiltonian describes a situation where the particles of A do not interact to each other, the same holds for the 
particles of B, but each particle of A interacts with all the particles of B and vice versa, as shown in Figure 6. 
In eq. (|41|). H is written in its diagonal form; then, the energy eigenvectors are 



i )...|^M-l)|^M)|Tl)...|tfc)-|TjV-l)|T J v) 

i> — 1-frAf-i) l-n-Ar) |Tl> — ]Tfc> — ITjv-i> Uiv) 



W-i) Wm) Ui) ... \lk) ... IliV-l) \In) 



(42) 



In turn, the eigenvectors of Ha form a basis of Ha- In order to simplify the expressions, we will introduce a particular 
arrangement into the set of those vectors, by calling them \Ai): the set {|-4i)} is an eigenbasis of Ha with 2 M elements. 
The \Ai) will be ordered in terms of the number I £ No of particles of A having spin |JJ.). Then, we have that: 



• / = corresponds to the unique state with all the particles with spin 

\Ar) = \t,t,-AA) 



H A \Ai) = y 



(43) 



• 1 = 1 corresponds to the M states with only one particle with spin |-||). Since the order of the eigenvectors with 
the same eigenvalue will be irrelevant for the computations, we will order these states in an arbitrary way: 



with j = 2,3, ...,M+1 



H A \Aj) 



M — 2 



\A 3 



(44) 
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(M-l)iU 



• 1 = 2 corresponds to the 
an arbitrary way: 



states with two particles with spin |-||). Again, we will order these states in 



\A 3 ) = \fa fa ,...AAA, , -AAA, -A A) Ha \Aj) = \Aj) 



with j -M + 2,M + 3, ...,M 



(M — 1) M 

1 + rr 1 



(45) 



• For the remaining values of I, the procedure is analogous. 



Consequently, we have: 



1 eigenvector with eigenvalue 



M 



M eigenvectors with eigenvalue 



M — 2 



Ml . . , . M — 21 
eigenvectors with eigenvalue 



(M — l)ll 



(46) 



with I = 0, 1, ...M. Then, it is clear that Ha is degenerate: it has 2 M eigenvectors but only M different eigenvalues. 
Therefore, a generic state \A) of the system A can be written in the basis {|-4j)} as 



\A)=Y,C i \A i )eH A with £|C,| 2 = 1 



i=l 



By introducing eq. (|4T|) into eq. (|4"0"|) , a pure initial state of the composite system U = A U B reads 



N 



■(a k \1k)+P k \lk)) 



vfc=l 



If we group the degrees of freedom of B in a single ket |£>(0)), |^ ) results 



(47) 



(48) 



\^ o )=Y,C l \A l )®\B(0)) 
i=i 



(49) 



The time-evolution of is ruled by the time-evolution operator U{t) = e lHt = e l ( H A®H B )t : 



\m) = «(t)iv»o> = E c * e-*^®^)* ia) ® ib(o)> - e c * e_lffA * 



-iH B t 



|B(0)> (50) 



If we use Afc to denote the eigenvalue of -Ha corresponding to the eigenvector |,Afc), then 



\m) = E ^ I- 4 *) ® e"^"*' |B(0)> = E °i I A> <8> |B(t)) 



(51) 



where (see eq. (|3T])) 



|B(t)) = e -' l ^ HBt \B(Q)) 



exp 



TV 



-iA k J29A\W(W-\lj)(li\)t 



|B(0)) 



(52) 
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Since the number of the eigenstates of Ha with the same eigenvalue is given by eqs. (|46)) . the terms of |^>(i)) can be 
arranged as 

/M+l \ /M+l+L^L 

\m) = (Ci IA) |fl (t))) + E C A |.Aa) |Bi(t)) + E C A |B 2 (i)) 

V A=l / V A=M+2 

+ E C ^ ^> + - + ( C 2M l<M*)» (53) 

where 

N 

Wi(t)) = (a*e <fli ^*« |tfc) + /V - *^*** |U>) (54) 
fc=i 

If we compare eq. (|53jl with eq. l[To]l . we can see that |£fr(t)) and are the particular cases of \Bi(t)) for M = 1 

and, then, I = 0, 1. Let us recall that Z is the number of particles of the system A having spin |JJ.). Then, with M = 1 
and I = 0, \Bi(t)) = \Sf(t)), and with M = 1 and I = 1, \Bi(t)) = \£^{t)). 
If we define the function 

then eq. ([55)) can be rewritten as 

M /(/) 

W*)>=E E Ca|^a)|^(*)) (56) 

i=0 A=/(i-l)+l 

and the state operator p(t) = \ijj{t)) (ip(t)\ reads 

/(!) 

M f(i') 

P(*)=E E C A Cjt,|^ A )|B,(*))(Bi'(*)IMv| (57) 



I,Z'=0 A = /(i-l) + l 
A' = /(('-l) + l 



B. Computing the expectation values 

An observable OEO = H<E)'Hoi the closed system U — A U B can be expressed as 

JV 



0= [ E S *,A< l- 4 ^ (^vl J ® ( (g) (4? IT,) (Til + eg |T,) (l,| + eg || 4 ) (Til + eg || 4 ) <| 4 |) j (58) 

Let us notice that eq. (|17p (a generic observable in the original spin-bath model) is a particular case of this eq. 
(|58p. with only four terms in the first factor. Analogously to that case, the diagonal components sa,a, eg, eg are 

real numbers, and the off-diagonal components are complex numbers satisfying sa.a' — s *\> a> e g = e g*' Then, the 
expectation value of the observable O in the state p(t) of eq. (|Ff|) can be computed as 

/(!) 

M /((') 

{0) p(t) =Tr(Op(t)) = E E B x,X'TiAt) (59) 

U'=0 A = /(!-l) + l 
A' = /(i'-l) + l 

where 

JV 



T I( j/(t) = J ! [lajfej^e^w-'-^'Oi + |^.| 2 e jJ e -*(9i.< + 2Re (a^^e't^+^'H )1 (60) 
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and 



g 3 .i = (21 - M) 9j , Byy = CxC* yS y. x (61) 

Since the exponents in cq. (|60[) are of the form g^i ± g^y , in some cases they are zero. So, we can write 



M 



f(M-l) 



M 



fW 
/('') 



M f(l) 

(0>K*)=E E ByyT u (t)+J2 E Byy2Re(T lM ^(t))+ Y, E ByyT^At) (62) 



1 = A = /(i-l) + l 
V=/(!-l) + l 



where 



1=0 A = /(!-l) + l 

A' = /(M-(-l) + l 



l,l'=0 A=/(l-l) + l 

!#(' a' = /(;-i) + i 



M 



M-2 



M-l 



if M is evenl 
if M is odd J 



(63) 



N 

Tun) = n [i 



a j\ e fT 



|/?/e«+2Re(a^e«e^ 



(64) 



Tl.M-l 



N 

(*) = II 4?^*" + l^l 8 4? e " i9,vt + 2Re (V 



.7=1 



(65) 



Let us notice that eqs. (|64|) and (|65|) are analogous to eqs. (|2Tj) and (|22]) for To(t) and Ii(t), respectively, in the 
original model, with g^i = (21 — M) gj instead of gj. In particular, when M = 1 and, so, Z = 0,1, then Tu(t) = To(t) 
andT,, M _,(t)=r 1 (t).' 

As in the case of the original spin-bath model, here we will consider different meaningful ways of selecting the 
relevant observables. 



VII. GENERALIZED SPIN-BATH MODEL: DECOMPOSITION 1 

A. Selecting the relevant observables 

In this case A is the open system S and B is the environment E. This is a generalization of Decomposition 1 in 
the original spin-bath model. The only difference with respect to that case is that here the system S is composed of 
M > 1 particles instead of only one. Then, the TPS for this case is 



Therefore, the relevant observables Or of the closed system U are those corresponding to A, and they are obtained 
from eq. (|55|) by making ei l ) = ef) = 1, ei l ) = (compare with eq. (f2"4"| in the original spin-bath model): 

Or = Os <8> 1 E = 

s\,\> \ A\) (Ay\ ® I (67) 

\A,A'=0 / \»=1 / 

With this condition, the expectation values of these observables are given by eq. (|6"2"|) . with 

N 

n 

N 



Ti,i(t) = n(i°ii a +i0ji j =i ( 6s ) 

N 

Ti, M -i(t) = n(l^| 2em!t + |/ 3 / e_m! *) ( 69 ) 



i=i 

N 



T U'(t) = U (|o! i | a e 4 («- , - ff J.")l + j^.^e^a*-'-^.'')*) (70) 
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If we define the functions Ri(t) = \T l:M -i{t)\ 2 and R w (t) = |T M /(i)| 2 , they result 

N 
N 

We can see that |r(i)| 2 of eq. (f27|) in the original model is the particular case of Ri(t) for M = 1. 





4 


H2|a,| 2 






4 







(71) 
(72) 



B. Computing the behavior of the relevant expectation values 

The expectation value given by eq. ([62j) has three terms, (0 R ) p{t) = + S (2) + £ (3) , which can be analyzed 
separately: 



From eq. ([68]) . the first term reads 



M /(*) M f(l) 

s(1) = E E *a,a' = E E c a q sv ,a^s (1) w (73) 

1=0 A,A'=/(i-l)+l Z=0 A,A'=/(i-l)+l 



It is clear that this first term does not evolve with time. 
• The time-dependence of the second term is given by Ti,M-i(t)'- 

/(<) 

M f(M-l) 



S(2) W=E E BA,A'2ReCZU_,(t)) (74) 



Z=0 A = /(i-l) + l 

V=/(M-I-1)+1 



Then, in order to obtain the limit of this term, we have to compute the limit of Ri(t) = \Ti.M-i(t)\ 2 of eq. (|71[) 



2 I 1 2 

As in the case of the original spin-bath model, here we take \aj\ and /3-J as random numbers in the closed 



interval [0, 1], such that \a 3 \ 2 + \f3 3 \ 2 = 1. Then 

max (\a 3 \ 4 + \P 3 \ 4 + 2\a 3 \ 2 \(3 j\ 2 cos (2(21- M) gjtj) =1 (75) 
mm (\a 3 \ 4 + A 4 + 2 \a 3 \ 2 \l3 J 2 cos (2(21- M)g 3 t)) = (2\a j \ 2 -l) (76) 



t 



Therefore, 



l a i| 4 + \Pj | 4 + 2 \ctj \ 2 \P 3 -\ 2 cos (2 (21 — M) g 3 t) is a random number which, if t ^ 0, fluctuates 
between 1 and ^2 |aj| 2 — lj . Again, when the environment has many particles (that is, when N — > oo), the 

2 I 1 2 2 I 1 2 

statistical value of the cases \a 3 \ — 1, \[3j\ — 1, \a 3 \ — and \/3A = tends to zero. In this situation, eq. 
(I7T|) for ii((t) is an infinite product of numbers belonging to the open interval (0, 1). As a consequence, when 
N -> oo, Ri(t) -> 0. 

• The time-dependence of the third term is given by Ti^i(t): 

tm 

s(3) W= E E B x,x'T h v{t) (77) 



I,!'=0 A = /((-l) + l 
A' = /(('-l) + l 



with the restrictions on / and (": Z ^ Z' and V ^ M — I. As in the second term, we have to compute the limit 
of Rw(t) = \T^i'(t)\ 2 of eq. ([72"]) and, on the basis of an analogous argument, the result is the same as above: 
when N — > oo, R u < (t) -» 0. 
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If we want now to evaluate the limit of (On) P (t) for t — > oo, we have to compute the limits of the second and the 
third terms (since the first term, as we have seen, is time- independent). Here we have to distinguish three cases: 
M < N, M > N and M ~ N. 

Case (a): M < N 

This case is similar to Decomposition 1 in the original spin-bath model, since in both cases M <C N: the only 
difference is that in the original model M = 1 whereas here M > 1. 

In fact, we have seen that Ti. M- i(t) is analogous to Ti(t) in the original model. Moreover, Tij>(t) has the same 
functional form as Fi(t). In paper [131 ] it is shown that Ti(t) approaches zero for t — > oo. This means that we can 
infer that Ti t M-i(t) and Tiji(t) also approach zero for t -> oo. On the other hand, the terms £( 2 )(i) and £ (3) (i) 
are sums of less than M terms involving Ti t M-i(t) and Ti^(t). As a consequence, since in this case M is a small 
number, the sum of a small number of terms approaching zero for t — > oo also approaches zero: lim^oo £( 2 )(t) = 
and limt^oo £( 3 )(t) = 0. Therefore, 



lim (0 R ) p{t) = lim [s«(t) + zW(t) + £( 3 )(t)l = £«(*) 



(78) 



In other words, 



M /(/) M /(0 

Hm (O fl ) p(t) = £ Y, B ^'=H E c aQ Sv ,a = (Ofl)p. 

Z=0 A,A'=/(i-l) + l i=0 A : A'=/(i-l) + l 



(79) 



where is the final diagonal state of ?7. This result can also be expressed in terms of the reduced density operator 
p A of the system A as (see eq. (O): 



lim (0 R ) p{t) = (0 R ) Pr = lim {0 A ) PA (t) = (Oa) Pa , 

L — ^OO ~L — ' OC 



(80) 



In the eigenbasis of the Hamiltonian Ha of A, the final reduced density operator p A ^ is expressed by a 2 M x 2 M 
matrix: 



Pa* 



(Pl=0 











... \ 





Pl=l 








... 








Pl=2 





... 











Pl=3 


... 


V o 











Pi =M ) 



(81) 



where p l=0 — \C\\ and each p l is a matrix of dimension n^fzjm\ X hftwi ■ This result might seem insufficient 
for decoherence because, since the p 2 are matrices, p A% seems to be non completely diagonal in the eigenbasis of 
the Hamiltonian H A - However, we have to recall that all the states \Ai) with same I are degenerate eigenvectors 
corresponding to the same eigenvalue of H A ; then, the basis that diagonalizcs p A ^ (i.e., that diagonalizcs all the 
matrices p ; ) is an eigenbasis of H A . Summing up, the system S = A of M particles in interaction with its environment 
E = B of N 3> M particles decoheres in the eigenbasis of p At , which is also an eigenbasis of H A . 

If we want to compute the time-behavior of (O r ) p ^, we have to consider that is a sum of terms of the form 

{B\ y |ay| + B\ y |/3J ), that is, terms of the expectation value coming from the diagonal part of pit) in the basis of 

the Hamiltonian H . Therefore, if there is decoherence, the sum E nd (t) = S (2) + S (3) , involving the terms of {0 R ) p{t) 
coming from the non-diagonal part of pit), has to approach zero for t — > oo. 

In order to show an example of the time-behavior of (0 R ) p M, numerical simulations for T, nd (t) have been performed, 
with the following features: 

(i) s A%A = 1 (see eq. ©). 

(ii) The initial condition for S = A is selected as (see eq. (fiTTO : 



1*4) 



Then, from (i) and (ii), B 



A, A' 



(see eq. (foTjl ). 




(82) 




0.0002 0.0004 0.0006 0.0008 0.001 

FIG. 7: Figure 7: Evolution of H nd {t) for N = 10 3 , and M = 1 (solid line) and M = 10 (dash line), with to = 10~ 3 s. 



(iii) | cti | 2 is generated by a random-number generator in the interval [0, 1], and I/3J 2 is obtained as |/3J 2 

(iv) gi — 400-ffz: as explained above, the coupling constant in typical models of spin interaction. 
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(v) As in the original model, the time- interval [0,to] was partitioned into intervals At = io/200, and the function 
£ nd (t) was computed at times tk — kAt, with k = 0,1, 200. 

(vi) N = 10 3 , and M = 1 and M = 10. 
Figure 7 shows the time-evolution of Y; nd (t). 

This result shows that, as expected, a small open system S — A of M particles decoheres in interaction with a large 
environment E = B of N ^S> M particles. 

Case (b): M > N 

In this case, where the open system S = A has much more particles than the environment E = B, the argument of 
Case (a) cannot be applied: since now Y,( 2 \t) and T,( 3 \t) are no longer sums over a small number of terms, the fact 
that each term approaches zero does not guarantee that the sums also approach zero. In particular, if N = 1, then 
(see eq. ([70)) ) 



T M ,(i) = | ai |V( ffl ' ! - ffl ' ! 'H + l^e-^.'-si,*')! 



(83) 



which clearly has no limit for t — > oo. Nevertheless, it might happen that, with TV high but M much higher than TV, 
each term of the sums approaches zero. So, in order to know the time behavior of (Or) p ^, numerical simulations 
for T, nd (t) have been performed, with the same features as in the previous case, with the exception of condition (vi), 
which was taken as: 



(vi) M = 10 3 , and TV = 10 and N = 100. 




0.002 0.004 0.006 0.008 

FIG. 8: Figure 8: Evolution of E nd (i) for M = 10 3 , and N = 10 (dash line) and N = 100 (dot line), with to = 10~ 3 s. 



Figure 8 shows the time-evolution of T, nd (t) in this case. 

This result is also what may be expected from the EID perspective: when the open system S = A of M particles 
is larger that the environment E = B of N <C M particles, S does not decohere. 

Case (c): M ~ N 

In this case, where the numbers of particles of the open system S = A and of the environment E = B do not differ 
in more than one order of magnitude, the time behavior of {Or) pit) cannot be inferred from the equations. Numerical 
simulations have been performed, with the same features as in Case (b), with the exception of condition (vi), which 
was taken as: 

(vi) N = 10 3 , and M = 10 2 and M = 10 3 . 

Figure 9 shows the time-evolution of £™ d (t). 

Again, this result is not surprising from the viewpoint of the EID approach: if the environment E = B of N 
particles is not large enough when compared with the open system S = A of M particles, S does not decohere. 

Let us notice that, for N = 10 3 , the system S = A with M = 10 2 does not decohere (Figure 9), whereas it does 
decohere with M = 10 (Figure 7). This shows that, in the case of this decomposition, M <§; N means that N is at 
least two orders of magnitude higher than M. 

Summarizing results 

Up to now, in this Decomposition 1 all the arguments were directed to know whether the system A of M particles 
decoheres or not in interaction with the system B of N particles. But, given the symmetry of the whole system, the 
same arguments can be used to decide whether the system B oi N particles decoheres or not in interaction with the 
system A of M particles, with analogous results: B decoheres only when M » N; if M <§C N or M ~ N, B does not 
decohere. Therefore, all the results obtained in this section can be summarized as follows: 



(i) If M <C N, A decoheres and B does not decohere. 



0.8 



0.6 



0.4 



0.2 








0.001 

FIG. 9: Figure 9: Evolution of E" d (t) for N = 10 3 , and M = 10 2 (dash line) and M = 10 3 (solid line), with to = 12.10" 4 s. 

(ii) If M ^> N , A does not decohere and B decoheres. 

(iii) If M ~ N, neither A nor B decohere. 

This general conclusion completely agrees with the usual reading of the EID approach: the decoherence of an open 
system S is the result of its interaction with a very large environment E. This happens when S — A in (i) and when 
S = B in (ii), but it does not happen in (iii). 



VIII. GENERALIZED SPIN-BATH MODEL: DECOMPOSITION 2 



A. Selecting the relevant observables 



In this case we decide to observe only one particle of the open system A. This amounts to splitting the closed 
system U into two new subsystems: the open system S is, say, the particle Am with ket IfMt, ff, ff, JJ-), and the 
environment is E = (u^f/j -1 ^) U B = (u*^ 1 ^) U (u^ =1 Bk). Let us notice that the Decomposition 2 of the original 
spin-bath model is a particular case of this one, for N = 1 (see eq. (|50|) . where N plays the role of the M of this 
case). The TPS for this case is 

H = H S ®H E = (h Am )® U§§ n A )j » (® Wb *JJ (84) 

Therefore, the relevant observables Or of the closed system U are those corresponding to the particle Am'- 

O r = O s ®1e=( J2 s ^-\a) (a'\ j ® ( ( (g) lA ® ((g) I* J J (85) 

\a,a'=ir,4 / V\i=l / \fc=l // 



20 



It is easy to see that the relevant observables selected in this Decomposition 2 form a subspace of the space of the 
relevant observables selected in Decomposition 1: cq. (|55|) can be obtained from eq. (|67[) by making s\ t \' = 1 for 
A = A' and s\ t y — for A ^ A' in all the terms of the sum except for the terms corresponding to the particle Am- 
In order to simplify expressions, in this case it is convenient to introduce a new arrangement for the eigenvectors 



of the Hamiltonian Ha, by calling them 
will be ordered by analogy with the binary numbers: 



Ai\: the set j -4«^> j is an eigenbasis of Ha with 2 M elements. The 



Ai) = 



fr,...,fr,fr,fr,fr>, A 2 ) = \t,t,-,1\,1\AA), As) = \tt -AAAA) , 



Ao 



A,-AAAA) 



At 



\1tA,-AAAA), 



(86) 



Ai ) with even i have the spin M in the state I-IJ-), and the 



According to this arrangement, the 
the spin M in the state Iff). So, the relevant observables of eq. (|85p can be rewritten in terms of the 



Ai ) with odd i have 



Ai ) as 



Or 



,A=1 



A 



2 A 



A 



2 A 



A 



2 A 



A 2 \-i 



A 2 \-i)(A 



^2A 



A 2 x-i) (A 2 \- 




(87) 



B. Computing the behavior of the relevant expectation values 

Here the expectation values of the relevant observables are given by eq. (|62"|) . with Ti t i>(t), Ti t i(t) and Tj af_j(i) 
given by eqs. (|6H)) . and (|6"5)) respectively, but now replacing B\ t \i with B\,y, 



M f{M-l) 



M 



/(!) 
/(!') 



M /(0 _ 
(°k)p(*)=E E ^,a'+E E S A , v 2Re(T z , M _ J (t))+ ]T ^ B x ^T l<v {t) (88) 



2=0 A=/(i-i) + i 

A' = /(l-l) + l 



=0 A = /(!-l) + l 

A'=f fM-l-ll+1 



i,l'=0 A=/(!-l)+l 
i^i' A'=/(l'-l) + l 



where the -B>, y can be written in the basi 



is { -4a)} 



as 



C\CyS^ft if A is an even number and A' = A 

CxCyS-ft^, if A is an even number and A' = A — 1 

B\,\' = < C\CyS^ if A is an odd number and A' = A + 1 

CxCyS]}.^ if A is an odd number and A' = A 
otherwise 



According to eq. (|89|) . B\ y ^ only when 



(89) 



A' = A or A' = A ± 1 
Since A = f(l - 1) + 1 and A' = /(/' - 1) + 1, relations imply that 

l' = l or I' = l±l 



(90) 



(91) 



The expectation value given by eq. (|88|) has again three terms, (O)p(t) = + E^ 2 ^ + Y,( 3 \ which can be analyzed 
separately: 



From eqs. (|89")1 and (|9U)l . the first term reads 

M /(/) 2 M - X 
S(1) =E E E (\C 2 x\ 2 S n + \C2X + i\ 2 S u )^^(t) 

1=0 A=/(/-l) + l A=0 

Analogously to eq. (|T5|) of Decomposition 1, this first term does not evolve with time. 



(92) 
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The time-dependence of the second term is given by Ti t M-i(t)- But with the restrictions of eqs. (|90| and (|91[) . 
£( 2 ) has only two terms: 



/CO 

M f(M-l) 

S(2) W = E E S A)V 2Re(T i , M -K*))= (93) 

Z = A = /(Z-1) + 1 

\'=/(M-!-l)+l 

= C/( Mjzi _i) + iCJ ( Af-i 1 - | i 2 (%f + %Jt) 2 Re ^ r^-i , m±i (t)^J (94) 

Then, in order to obtain the limit of this term, we have to compute the limit of T m-i m+i (t), which is precisely 

the Ti t ii(t) of Decomposition 1 in the particular case that I — M ~ 1 and V — M ^ (see eq. ([70]) ). But, as we have 
seen in Case (a) of Decomposition 1, Tj has the same functional form as T\{t) of the original model (see eq. 
([2"2"| ). which approaches zero for t — ► oo when JV 3> 1. Therefore, for AT 3> 1, Tm-i m+i (t) also approaches zero 

for £ — > oo, and the same holds for T,^'(t) since it is a sum of two terms containing Tm-i m+i (t). 

• The time-dependence of the third term is given by Ty (t). But with the restrictions of eqs. (|90[) and (|9T]) . £( 3 ) 
results: 

M /(0 

S (3) W= E E (£A,A+iT M+ i(t) + 5 A ,A-iT M _i(i)) (95) 
Since here /' = I ± 1 (see eq. (|9"Tj)). in this case Tij±i(t) is: 

t m ±i(<) = n (i«ii 2e=Fi9it + i^r^*) (96) 

If we compare this equation with eq. (|26p for r(£) in the original spin-bath model, we can see that 

Ti,i +1 (t) = r(t) and T M _i(t) = r*(t) (97) 

Then, 

E( 3 Hf) = (5 + r(t) + S_r*(t)) (98) 

where and SL are constants given by 

A/ /(/) 

s±= E E B ^±i (") 

'5°-! *=/0-i)+i 

,?: 2 

On the basis of the simulations of the original model we have seen that, when N ^S> 1, r(t) approaches zero for 
t — ► oo. Therefore, in this case we can conclude that, when JV 3> 1, £( 3 )(J) approaches zero for £ — > oo. 

Summing up, {On) P (t) is the sum of three terms: one is time-independent and the other two tend to zero for t — > oo. 
In particular, from eq. (j9"2"|) we know that, for JV 3> 1, 

M /(/) m 2 M ~ 1 

hm (O fl ) p(t) =^ £ fl A , v = E E (I^aI 2 % + |C 2 a+i| 2 S M ) = <O fl ) P . (100) 

Z=0 A,A'=/(/-l)+l '=0 A=0 

where p + is the final diagonal state of £7. Again, this result can also be expressed in terms of the reduced density 
operator p s = p Au of the open system S = Am as (see eq. ([80])) 

hm (0 R ) p(t) = (0 R ) Pr = hm (0 Am ) p . (t) = (Ox*),. . (101) 
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FIG. 10: Figure 10: Evolution of E nd (i) for M = 10 3 and N = 1, with to = 3.10" 2 s. 



where the final reduced density operator p Am * in the basis {Iff) , |JJ-)} reads 

This shows that the open system S = Am, composed of a single particle, decoheres in interaction with its environment 
E of N + M — 1 particles when N ^> 1, independently of the value of M. 

In order to illustrate this conclusion, we have computed Yi nd {t) = (t) + Y,^ (t) by means of numerical simulations 
with the same features as in Decomposition 1, with the exception of condition (vi), which was taken as: 

Figure 10: (vi) M = 10 3 and N = 1. 

Figure 11: (vi) M = 10 3 and N = 10 2 . 

Figure 12: (vi) M = 10 3 and N = 10 3 . 

Summarizing results 

As we have seen, in this decomposition of the whole closed system, the open system S = Am decoheres when 
N ^> 1, independently of the value of M. But the particle Am was selected as S only for computation simplicity: the 
same argument can be developed for any particle Ai of A. Then, when N ^> 1 and independently of the value of M, 
any particle Ai decoheres in interaction with its environment E of N + M — 1 particles. 

On the other hand, as in Decomposition 1, here the symmetry of the whole system U allows us to draw analogous 
conclusions when the system S is one of the particles of B, say, B^: S = Bn decoheres when M ^> 1, independently 
of the value of N. And, on the basis of the same considerations as above, when M 3> 1 and independently of the 
value of N, any particle B t decoheres in interaction with its environment E of N + M — 1 particles. 




IX. DISCUSSION 
A. Analyzing results 

According to the usual reading of the EID approach, the decoherence of an open system is induced by its interaction 
with a large environment. Such an interaction is what leads to the dissipation of energy from the open system S to the 
environment E. So, the orthodox view suggests a picture of decoherence where an energy flow from the open system 
S to the environment E washes out the original coherence and allows the classicality of S to emerge. According to 
this picture, in the original spin-bath model a spin-1/2 particle S = P decoheres when immersed in a large bath of 
spin-1/2 particles E = UiPi: P dissipates its energy into E and may acquire a classical nature. However, this reading 
has to face the "looming big" problem of defining the open systems involved in decoherence, since it does not provide 
a criterion to identify the open system S and its environment E. Now we will discuss the results obtained in the 
generalized spin-bath model, in order to see how they may contribute to the clarification of the problem. 

a) As we have seen, in our generalized model, where U = A U B, with A of M particles Ai and B of N particles 
Bi, (i) when M > JV or M ~ JV, the subsystem A does not decohere (Decomposition f of Section VII), but (ii) the 
particles A i} considered independently, decohere when N » 1 (Decomposition 2 of Section VIII). This means that 
there are physically meaningful situations, given byM>JV>lorM~JV>l, where all the Ai decohere although 
A does not decohere. In other words, in spite of the fact that certain particles decohere and may behave classically, 
the subsystem composed by all of them retains its quantum nature. This seemingly paradoxical conclusion sounds 
even more strange when the situation is conceived in terms of energy dissipation. In spite of the fact that all the Ai 
dissipate their energy into the environment (mainly into the subsystem B due to the interaction among each A t and 
all the Bi), the composite system A = UiAi (which should dissipate the energy of all the Ai) does not decohere. 

b) We have also seen that, by symmetry, all the particles Bi, considered independently, also decohere when M ^> 1. 
Then, when M>JV>lorM~JY>l, the requirement M ^> I holds and we can conclude that not only all the 
Ai, but also all the Bi decohere. So, all the particles of the closed system U = (Uj-Aj) U (DjBj) may become classical 
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FIG. 12: Figure 12: Evolution of E nd (i) for M = 10 3 and N = 10 3 , with to = 4.1(T 4 s. 

when considered independently, although the whole system U certainly does not decohere and, therefore, retains its 
quantum character. Again, the explanation of this result is even more difficult when it is conceived in terms of the 
energy dissipated from the system that decoheres to its environment, since we are committed to decide which particles 
give and which receive the dissipated energy. 

These difficulties are further consequences of the "looming big" problem of defining the open systems involved in 
decoherence. The irony of this story is that such a problem is the consequence of what has been considered to be the 
main advantage of the decoherence program: its open-system perspective. According to this perspective, particles 
interacting with other particles by exchanging energy are well-defined open systems, and the collections of those 
particles are open systems too. So, the problem is to decide which one of all these open systems is the system that 
decoheres or, in other words, where to place the cut between the system S and its environment E. 

The open-system approach not only leads to the "looming big" problem, but also disregards the well-known holism 
of quantum mechanics: a quantum system in not the mere collection of its parts and its interactions. In order to 
retain its holistic nature, a quantum system has to be considered as a whole: the open "subsystems" are only partial 
descriptions of the whole closed system, given by the selection of particular subspaces of relevant observables. On the 
basis of this closed-system perspective, we can develop a different conceptual viewpoint for understanding decoherence, 
which dissolves the problems of the orthodox open-system view. 

B. A different conceptual viewpoint 

As we have seen, a TPS expresses the decomposition of the closed system U, represented in the Hilbert space 
H = Ha <8> Hb, into two open systems Sa and Sb, represented in Ha and Hb respectively. Such a decomposition 
amounts to the split of the whole space O = H <8> H of the observables of U into the subspaces Oa = Ha <8> Ha and 
Ob = Hb <S> Hb such that O = Oa <8> Ob- In particular, the total Hamiltonian of U, H S 0, can be expressed as 
H = Ha ® Ib + 1 A® Hb + Hab , where Ha & Oa is the Hamiltonian of Sa , Hb € Ob is the Hamiltonian of Sb , and 
Hab € O is the interaction Hamiltonian, representing the interaction between the systems Sa and Sb- 
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As stressed in papers fTTJ] and [11], in general a quantum system U admits a variety of TPSs, that is, a variety 
of different decompositions into S A and Sb , each one defined by the space of observables Oa of Sa and Ob of Sb- 
Among all the possible decompositions of U, there is a particular TPS that remains dynamically invariant. This is 
the case when the interaction Hamiltonian Hab = 0: there is no interaction between Sa and Sb and, then, 

[H A ®I b ,Ia® H b ] = => cxp (-iHt) = exp (-iH A t) cxp (-iH B t) (103) 

Therefore, 

p A (t) = Tr (B) p(t) = e iH ^ (Tr (B)Po ) = e^Vxo e""™ (104) 

p fl (t) = Tr [A) p(t) = e iHBt (Tr {A) p ) **** = e iHBt p B0 e^ Bt (105) 

This means that, even if the initial state p of U is an entangled state with respect to the TPS H — Ha ® Hb, the 
subsystems Sa and 5s are dynamically independent: each one of them evolves unitarily under the action of its own 
Hamiltonian. As a consequence, the subsystems S A and Sb resulting from this particular, dynamically invariant TPS 
do not decohere. 

Once we have excluded the dynamically invariant TPS of U , all the remaining TPSs define interacting subsystems 
S A and Sb, such that Hab 7^ 0. As a result of the interaction, Sa and Sb evolve non-unitarily and, then, depending 
on the particular interaction between them, they may decohere. But the point to stress here is that there is no 
privileged non- dynamically invariant decomposition of U: each partition of the closed system into Sa and Sb is just 
a way of selecting the spaces of observables Oa and O b ■ 

When we adopt a closed-system perspective by means of the concept of TPS, it turns out to be clear that, in 
decoherence, there is no essential criterion for identifying the "open system" and its "environment" . Given the closed 
system U, that identification requires two steps: (i) to select a TPS H = Ha®Hb, such that U — Sa^Sb, and (ii) to 
decide that one of the systems resulting from the decomposition, say S A , is the open system S, and the other, S B , is the 
environment E. Since the TPS is defined by the spaces of observables Oa and Ob, the decomposition of U is just the 
adoption of a descriptive perspective: the identification of S and E amounts to the selection of the observables relevant 
in each situation. But since the split can be performed in many ways, with no privileged or essential decomposition, 
there is no need of an unequivocal criterion for deciding where to place the cut between "the" system and "the" 
environment. Decoherence is not a yes-or-not process, but a phenomenon relative to the chosen decomposition of the 
whole closed quantum system. When viewed from this closed-system perspective, Zurek's "looming big problem" does 
not constitute a real threat to the decoherence program: the supposed challenge dissolves once the relative nature of 
decoherence is taken into account. 

^,From this perspective, quantum mechanics is a theory whose dynamical postulate refers to closed systems: the 
time-behavior of the parts resulting from different partitions of the closed system has to be inferred from that postulate. 
Since the total Hamiltonian rules the dynamical evolution of the closed system, then the time-behavior of its open 
subsystems depends on the form in which the Hamiltonian is decomposed in each particular partition. This means that 
decoherence cannot be simply described as the result of an interaction through which a small open system —typically, 
a particle— dissipates its energy into a large environmental bath. As we have seen in the generalized spin-bath model, 
this picture of decoherence leads to perplexities: the relationships between the whole closed system and its open 
subsystems is subtler than that picture suggests (in a future paper we will study those relationships from a theoretical 
viewpoint in order to draw some general conclusions regarding decoherence). Therefore, the decomposition of the 
total Hamiltonian has to be studied in detail in each particular case, in order to know whether the system of interest 
resulting from the partition of the whole closed system decoheres or not under the action of its self-Hamiltonian and 
the interaction Hamiltonian. 



X. CONCLUSIONS 

The aim of this paper has been to argue that decoherence can be viewed from a closed-system perspective, which 
improves the understanding of the phenomenon. For this purpose, we have analyzed the simple spin-bath model by 
studying the time-behavior of the expectation values of relevant observables belonging to different sets. Then, we have 
generalized the original model in order to see how decoherence depends on the way in which the relevant observables 
are selected. 

On the basis of the analysis of the two models from a closed-system perspective, we have drawn the following 
conclusions: 

(i) Decoherence is a phenomenon relative to which degrees of freedom of the whole closed system are considered 
relevant and which are disregarded in each case. 
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(ii) The explanation of decoherence requires the detailed study of the interaction Hamiltonian resulting from the 
selected partition of the whole closed system. 

(ii) Although it is usually claimed that EID is a dissipative approach to decoherence, the simple account of deco- 
herence in terms of energy dissipation from the open system to its environment is misguided, to the extent that 
there are situations where all the particles of a closed system decohere when considered independently. 

(iv) Since there is no privileged or essential decomposition of the closed system, there is no need of an unequivocal 
criterion for identifying the systems involved in decoherence. Therefore, the "looming big problem" , which, 
according to Zurek, poses a serious threat to the whole decoherence program, looses its strength in the light of 
the relative nature of decoherence. 
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